A Parallel Unsymmetric Inverse Iteration Solver
نویسنده
چکیده
We describe a matrix multiply based block unsymmetric inverse iteration solver for upper Hessenberg matrices. Our kernel is robust in that it prevents overrow by scaling. It uses new techniques to ensure performance is not sacriiced when scaling is not necessary. Finally, we give results on a parallel implementation on an Intel Paragon TM supercomputer.
منابع مشابه
An Approach for Parallelizing any General Unsymmetric Sparse Matrix Algorithm
In many large scale scientiic and engineering computations, the solution to a sparse linear system is required. We present a partial unsymmetric nested dissection method that can be used to parallelize any general unsymmetric sparse matrix algorithm whose pivot search can be restricted to a subset of rows and columns in the active submatrix. The submatrix is determined by the partial unsymmetri...
متن کاملA CPU-GPU hybrid approach for the unsymmetric multifrontal method
Multifrontal is an efficient direct method for solving large-scale sparse and unsymmetric linear systems. The method transforms a large sparse matrix factorization process into a sequence of factorizations involving smaller dense frontal matrices. Some of these dense operations can be accelerated by using a graphic processing unit (GPU). We analyze the unsymmetricmultifrontalmethod fromboth an ...
متن کاملNUMERICAL ANALYSIS GROUP PROGRESS REPORT January 1994 – December 1995
2 Sparse Matrices ……………………………………………………………………………… 4 2.1 The direct solution of sparse unsymmetric linear sets of equations (I.S. Duff and J.K. Reid) …………………………………………………………………………… 4 2.2 The design and use of algorithms for permuting large entries to the diagonal 2.6 Element resequencing for use with a multiple front solver (J. A. Scott) ………… 10 2.7 Exploiting zeros on the diagonal in the direct s...
متن کاملUnsymmetric-pattern Multifrontal Methods for Parallel Sparse Lu Factorization
Sparse matrix factorization algorithms are typically characterized by irregular memory access patterns that limit their performance on parallel-vector supercomputers. For symmetric problems, methods such as the multifrontal method replace irregular operations with dense matrix kernels. However, no e cient method based primarily on dense matrix kernels exists for matrices whose pattern is very u...
متن کاملIterative Methods for Neutron Transport Eigenvalue Problems
We discuss iterative methods for computing criticality in nuclear reactors. In general this requires the solution of a generalised eigenvalue problem for an unsymmetric integro-differential operator in 6 independent variables, modelling transport, scattering and fission, where the dependent variable is the neutron angular flux. In engineering practice this problem is often solved iteratively, u...
متن کامل